3. Greens functions. H = H 0 + V t. V t = ˆBF t. e βh. ρ 0 ρ t. 3.1 Introduction
|
|
- Neal Park
- 5 years ago
- Views:
Transcription
1 3. Greens functions 3.1 Introduction Greens functions appear naturally as response functions, i.e. as answers to the function how a quantum mechanical system responds to an external perturbation, like for example electrical or magnetic fields; the corresponding response functions would then describe the electrical conductivity or the magnetic susceptibility of a system. Here we will be concerned with small perturbations and thus only the linear response of the system. We describe the system by a Hamiltonian H = H 0 + V t (3.1) where V t represents the interaction with an external field. H 0 describes the system with the external field switched off; due to interactions H 0 is not necessarily exactly solvable. The external field F t couples to the observable ˆB of the system: V t = ˆBF t (3.2) Here, ˆB is an operator and F t is a complex number. We now consider an observable  of the system that is not explicitly time dependent and ask how the dynamic expectation value  reacts to the perturbation V t. Without field we have  0 = Tr(ρ 0Â) (3.3) where ρ 0 is the density matrix of the system without external fields: ρ 0 = e βh Tre βh (3.4) in the grand canonical ensemble H = H µ ˆN (with chemical potential µ, particle number operator ˆN). The density matrix will change if we switch on the field: ρ 0 ρ t (3.5) 29
2 This means for the expectation value of   t = Tr(ρ tâ) (3.6) In the Schrödinger picture, the equation of motion of the density matrix (the statistical operator) is determined by the von Neumann equation i h ρ t t = [ H + V t, ρ t ] (3.7) We consider a perturbation that is switched on at some time so that the boundary condition for our first order differential equation is an unperturbed system for t lim ρ t = ρ 0. (3.8) t We now switch to the Dirac picture where we have ρ D t (t) = e i h H 0t ρ t e i h H 0t (3.9) with the equation of motion ρ D t (t) = i h[ ρ D t, V D t Integrating with the boundary condition leads to ] (t). (3.10) lim t ρd t (t) = ρ 0 (3.11) ρ D t (t) = ρ 0 i h t dt [ V D t (t ), ρ D t (t ) ] (3.12) This equation can be solved by iteration (by substituting ρ D t (t) repeatedly on the right hand side): ρ D t (t) = ρ 0 + ρ D(n) t (t) = n=1 ( i h ) n t ρ D(n) t (t) with dt 1 t1 tn 1 dt 2 dt n [ V D t 1 (t 1 ), [ V D t 2 (t 2 ), [ [V D t n (t n ), ρ 0 ] ] 30 ] ] (3.13)
3 While this formula is exact, it is not practical. For sufficiently small external perturbations, we can restrict to linear terms in the perturbation V t which is called linear response: ρ t ρ 0 i h t dt e i h H 0t [ V D t (t ), ρ 0 ] e i h H 0t (3.14) Here, we have returned to the Schrödinger representation for the density matrix. We can use this result to determine the perturbed expectation value of (3.6):  t = Tr(ρ tâ) =  0 i h =  0 i h =  0 i h t t t dt Tr {e [ ] i h H 0t Vt D (t ), ρ 0 { ] dt F t Tr [ˆB D (t ), ρ 0 } e i h H 0t ÂD (t) =ˆBρ0  ρ 0ˆBÂ=ρ 0 ˆB ρ0ˆbâ { dt F t Tr ρ [ÂD 0 (t), ˆB D (t ) ] } } (3.15) Here, cyclic invariance of the trace was exploited. This shows how the system reacts to the external perturbation, as measured from the observable Â: A t =  t  0 = i h t dt F t [ÂD (t), ˆB D (t )] 0 (3.16) This response is determined by an expectation value of the unperturbed system. The Dirac representation of the operators ÂD (t), ˆB D (t ) corresponds to the Heisenberg representation when the field is switched off. Now we define the retarded two-time greens function G r AB(t, t ) = A(t); B(t ) = iθ(t, t ) [ A(t), B(t ) ] (3.17) 0 The operators are to be taken in Heisenberg representation of the field free system. The retarded Greens function describes the response of a system as manifested in observable  when the perturbation couples to observable ˆB: A t = 1 h dt F t G r AB(t, t ) (3.18) 31
4 It is called retarded because due to the Heaviside function, only perturbations for t < t contribute. With the Fourier transform F(ω) of the perturbation F t = 1 2π dω F(ω)e i(ω+iδ)t (3.19) where δ > 0 is infinitesimally small and using the later result that with a Hamiltonian that is not explicitly time dependent the Greens function depends only on time differences t t, we can rewrite (3.18) in the form of the Kubo formula A t = 1 dω F(ω)G r 2π h AB(ω + iδ)e i(ω+iδ)t (3.20) The δ > 0 in the exponent enforces the boundary condition (3.8). We will now look into two applications of response functions. Magnetic Susceptibility The perturbation is a spatially homogeneous magnetic field that oscillates in time: B t = 1 dω B(ω)e i(ω+iδ)t, (3.21) 2π which couples to the magnetic moment m = m i = gµ B S i. (3.22) h i i Thus, the perturbing potential term in the Hamiltonian becomes V t = m B t = 1 dω m α B α (ω)e i(ω+iδ)t (3.23) 2π α where α = x, y, z are Cartesian directions. An interesting quantity is now the magnetization in response to the applied field. As it is M = 1 V m = gµ B S i, (3.24) hv i we have to choose the magnetic momentum operator for both  and ˆB operators in the Kubo formula: M β t Mβ 0 = 1 dt B α t m β (t); m α (t ). (3.25) V α 32
5 Only in a ferromagnet there is a finite magnetization M β 0 Eq. (3.25) defines the magnetic susceptibility tensor without a field. χ βα ij (t, t ) = µ 0 gµ 2 B S β V h 2 i (t); Sα j (t ) (3.26) as a retarded Greens function. Thus M β t = 1 dt B α t µ χβα ij (t, t ) (3.27) 0 ij or in terms of frequency M β t = 1 2πµ 0 ij α dω e i(ω+iδ)t χ βα ij (ω)b α (ω) (3.28) We have implicitly assumed that the system we consider has permanent localized moments. Two types of susceptibilities are interesting: The longitudinal susceptibility χ zz ij (ω) = µ 0 gµ 2 B S z V h 2 i ; S z j (3.29) ω where the index indicates the Fourier transform of the retarded Greens function. This can be used to obtain information about the stability of magnetic orderings. For the paramagnetic phase, one calculates the spatial Fourier transform χ zz q (ω) = 1 N ij ( Ri ) χ zz ij (ω)e i q R j (3.30) At the singularities of this response function, an infinitesimally small field is sufficient to create a finite magnetization, i.e. a spontaneous ordering of the moments. For that purpose, the conditions under which { lim ( q,ω) χ zz q (ω)} 1 = 0 (3.31) are studied; they indicate the paramagnetic ferromagnetic transition. The other interesting case is the transversal susceptibility χ + ij (ω) = µ 0 gµ 2 B S + V h 2 i ; S j ω where S± i = S x i ± is y i (3.32) 33
6 Poles of this susceptibility correspond to spin wave (magnon) energies: { } 1 (ω) = 0 ω = ω( q). (3.33) χ + q The examples show that linear response theory not only treats weak external perturbations but also yields information about the unperturbed system. Electrical conductivity Now we consider a spatially homogeneous electrical field that oscillates in time: E t = 1 2π dω E(ω)e i(ω+iδ)t. (3.34) The electrical field couples to the electrical dipole moment P P = d 3 r rn( r). (3.35) We consider N point changes q i at positions r i (t); the charge density is n( r) = N q i δ( r r i ). (3.36) i=1 This gives a dipole moment operator P = N i=1 q i r i. (3.37) The electrical field causes the additional external potential term in the Hamiltonian V t = P E t = 1 dω P α E α (ω)e i(ω+iδ)t. (3.38) 2π α An interesting quantity is the response of the current density to the external field: j = 1 V N i=1 q i r i = 1 P. (3.39) V 34
7 Its expectation value without field disappears: j 0 = 0. (3.40) After switching the field on, we have j β t = 1 h dt E α t j β (t); P α (t ). (3.41) α In terms of the Fourier transforms this becomes j β t = 1 dω e i(ω+iδ)t σ βα (ω)e α (ω) (3.42) 2π α This is Ohms law, defining the electrical conductivity tensor σ βα (ω) j β ; P α ω (3.43) that has retarded Greens functions as components. This can be rewritten as σ βα (ω) = i N q 2 V m(ω + iδ) δ αβ + i h j β ; j α (3.44) ω + iδ The first term represents the conductivity of a noninteracting electron system as given by classical Drude theory, and the second one involving a retarded current-current Greens function represents the interaction between the particles. 3.2 Matsubara method In the solid state theory class, Greens functions were introduced as response functions; they can be used to determine the quasiparticle density of states. They also appear as correlation functions, and they give us access to excitation energies. But so far, everything was done at zero temperature. The ojective is now to extend the methods for T = 0 to finite temperatures. 1 As a reminder, the two-time Greens functions were defined as (with ε = 1 for fermionic operators, ε = +1 for bosonic operators) retarded: G r AB(t, t ) A(t); B(t ) r = iθ(t t ) [A(t), B(t )] ε advanced: G a AB(t, t ) A(t); B(t ) a = iθ(t t ) [A(t), B(t )] ε causal: G c AB(t, t ) A(t); B(t ) c = i T ε (A(t)B(t )) 1 This chapter is based on K. Elk, W. Gasser, Die Methode der Greenschen Funktionen in der Festkörperphysik, Akademie-Verlag Berlin 1979, and W. Nolting, Grundkurs Theoretische Physik 7, Viel-Teilchen-Theorie, Springer
8 with the Wick time ordering operator T ε : (3.45) T ε (A(t)B(t )) = θ(t t )A(t)B(t ) + εθ(t t)b(t )A(t) (3.46)... indicates an average over the grand canonical ensemble X = Tr(e βh X) Tr(e βh ) where β = 1, H = H µn (3.47) k B T If the Hamiltonian has no explicit time dependence, the Greens functions are homogeneous in time: H t = 0 Gα AB(t, t ) = G α AB(t t ), α = r, a, c (3.48) This can be shown using the cyclic invariance of the trace. For the correlation functions that are needed for G r and G a, we have (A(t)B(t ) = A(t t )B(0) = A(0)B(t t) (B(t )A(t) = B(t t)a(0) = B(0)A(t t ) (3.49) We now allow the time variables to formally take complex values: Tre βh A(t i hβ)b(t ) = Tr { e βh e i h H(t i hβ) A(0)e i h H(t i hβ) B(t ) } = Tr { B(t )e βh e +βh e i h Ht A(0)e i h Ht e βh} = Tr { e βh B(t )A(t) } A(t i hβ)b(t ) = B(t )A(t) (3.50) As two different correlation functions become related in this way, the extension of the Greens function to complex time seems to be advantageous. In particular, in perturbation theory in V where H = H 0 + V, V would, for finite temperatures, appear in two places, in the Heisenberg representation of time dependent operators e ± i h Ht and in the density operator of the grand canonical averaging e βh ; two perturbation expansions would be necessary. Therefore, we join the exponential functions by introducing a complex time. The Matsubara method introduces purely imaginary times so that the quantity τ = it is real. This leads to a modified Heisenberg representation of operators: A(τ) = e 1 h Hτ A(0)e 1 h Hτ (3.51) 36
9 Note that the operator e i h Hτ creating imaginary time shifts is not unitary. The equation of motion (EOM) for an operator A(τ) h τ A(τ) = h τ [e 1 h Hτ A(0)e 1 h Hτ ] = HA(τ) + A(τ)H (3.52) thus becomes: h τ A(τ) = [ A(τ), H ] (3.53) Here, we use the (conventional) step function θ(τ) = { 1 for τ > 0 (t = τi negative imaginary) 0 for τ < 0 (t = τ i positive imaginary) (3.54) It can be used to introduce the time ordering operator. T τ { A(τ)B(τ ) } = θ(τ τ )A(τ)B(τ ) + ε p θ(τ τ)b(τ )A(τ) (3.55) where p is the number of permutations of creation operators. We assume pure Fermi/Bose operators so that p = 1. The definition of the Matsubara function is: G M AB(τ, τ ) = A(τ); B(τ ) M = T τ (A(τ)B(τ )) (3.56) Using (3.53) and (3.55) we get the EOM h τ GM AB(τ, τ ) = hδ(τ τ ) [A, B] ε + [A(τ), H] ; B(τ ) (3.57) Properties of the Matsubara Greens function introduced in this way are: 1) it depends only on time differences G M AB(τ, τ ) = G M AB(τ τ, 0) = G M AB(0, τ τ) (3.58) 37
10 2) it is periodic in time: Take hβ > τ τ + n hβ > 0, n Z; then Tre βh G M AB(τ τ + n hβ) >0 = Tr { e βh T τ ( A(τ τ + n hβ)b(0) )} = Tr { e βh A(τ τ + n hβ)b(0) } = Tr { e βh e H h (τ τ +n hβ) A(0)e H h (τ τ +n hβ) B(0) } = Tr { e H h (τ τ +(n 1) hβ) A(0)e H h (τ τ +(n 1) hβ) e βh B(0) } = Tr { } e βh B(0)A(τ τ + (n 1) hβ) = Tr { e βh T τ ( B(0)A(τ τ + (n 1) hβ) )} <0 = εtr { e βh T τ ( A(τ τ + (n 1) hβ)b(0) )} (3.59) This gives us the important result: ( hβ > τ τ +n hβ > 0 : τ τ +n hβ ) = εgab( M τ τ +(n 1) hβ ) In particular, for n = 1: G M AB (3.60) G M AB( τ τ + hβ ) = εg M AB( τ τ ) for hβ < τ τ < 0 (3.61) Thus the Matsubara Greens function is periodic in an interval 2 hβ; it is enough to consider the time interval hβ < τ τ < 0. This periodicity allows the Fourier expansion: G M (τ) = 1 [ 2 a 0 + a n cos nπ hβ τ + b n sin nπ ] hβ τ a n = 1 hβ b n = 1 hβ n=1 hβ hβ Now we define discrete energies E n = nπ β ( nπ ) dτ G M (τ) cos hβ τ dτ G M (τ) sin ( nπ hβ τ ) (3.62) (3.63) and the Matsubara Greens function on the imaginary energy (frequency) axis G(E n ) = 1 2 hβ(a n + ib n ) (3.64) 38
11 Then G M (τ) = 1 hβ G M (E n ) = 1 2 e i h Enτ G M (E n ) n= hβ dτ G M (τ)e i h E nτ (3.65) Still, we can simplify a bit more: G M (E n ) = 1 2 = dτ G M (τ)e i h E nτ dτ G M (τ)e i h E nτ = [ 1 + εe ie nβ ] hβ 0 dτ G M (τ)e i h E nτ dτ G M (τ)e i h E nτ dτ G M (τ hβ)e i h E nτ e ie nβ (3.66) where τ = τ + hβ was introduced. The bracket [... ] disappears for fermions (ε = 1) if n is even, for bosons (ε = 1) if n is odd. Thus with G M (τ) = 1 hβ G M (E n ) = E n = n= 0 e i h E nτ G M (E n ) dτ G M (τ)e i h E nτ { 2n π β for bosons (2n + 1) π β for fermions (3.67) (3.68) Now we need to work out a spectral representation for the Matsubara Greens function in order to relate it to the retarded Greens function. We first consider the correlation function 1 A(τ)B(0) = Tre βh E n A(τ)B(0) E n e βe n n 1 = Tre βh E n A E m E m B E n e βe n e 1 h (E n E m )τ (3.69) n 39
12 where we introduced the eigenstates of the Hamiltonian E n : H E n = E n E n. We recall the spectral density on the real energy axis: S AB (E) = h Tre βh E n A E m E m B E n e βe n (1 εe βe )δ ( E (E m E n ) ) n,m A(τ)B(0) = 1 h de S AB(E) 1 εe βee 1 h Eτ (3.70) In the integration interval in (3.67), τ is positive; we need to evaluate G M AB(E n ) = Now we use the integral 0 0 dτ e 1 h (ie n E)τ = dτ e i h E nτ A(τ)B(0) (3.71) h ( e iβe n e βe 1 ) = ie n E and put Eq. (3.70) into Eq. (3.71) to obtain G M AB(E n ) = h ( εe βe 1 ) ie n E (3.72) de S AB(E ) ie n E (3.73) If we compare this to the spectral representation of the retarded Greens function, we only need to replace ie n E+i0 + ; thus, the retarded Greens function can be obtained from the Matsubara Greens function by analytic continuation! 3.3 Some methods for Matsubara axis functions The Green s function on the Matsubara axis is fairly smooth and featureless so that the determination of the sums poses no problems. But there are two points to be considered: (1) Terms falling off as 1/(iω n ) will be badly represented by any number of frequency points; therefore, high frequency corrections are needed. (2) In order to obtain a spectral function, we need to analytically continue the Greens function to the real axis. This is done by the Padé method. 40
13 High frequency correction The idea is to subtract the leading terms of the high frequency expansion and calculate them analytically, e.g. F m (iω n ) = 1 n f + n f ( N 1 ) + O iω n (iω n ) 2 (3.74) n f = T (F m (iω n) 1 n f + n f N )e iω n0 + iω + 1 n f + n f N iω n n 2 (3.75) The last term of Eq. (3.75) is obtained by explicit calculation of the sum: 1 T e iωnτ 2i sin(ω n τ) = T iω n= n iω n=0 n = 2 sin [ (2n + 1)πTτ ] = 1 (3.76) π 2n n=0 because τ = 0 + ]0; 1/T[ and 4 sin [ (2n + 1)x ] 1 for π < x < 0 = 0 for x = π, 0, π π 2n + 1 n=0 1 for 0 < x < π In practice, the 1/(iω n ) 2 is taken into account in the same way. (2) Pade approximation This method for analytic continuation can be used (3.77) 1. when the function to be continued is not given analytically (in that case, use iω n ω + iδ) 2. when the function is given without statistical errors (for Quantum Monte Carlo Green s functions, use Maximum Entropy Method) One has to keep in mind one drawback of the Padé approximation: it is a polynomial representation that has limited precision for functions which are hard to approximate by a polynomial. It can show the wiggles that are typical for polynomial interpolations. The algorithm for calculating the Padé approximant was written down nicely by Vidberg and Serene (H. J. Vidberg and J. W. Serene, Solving the 41
14 Eliashberg Equations by Means of N-Point Padé Approximants, J. Low Temp. Phys. 29, 179 (1977)). Given the values u i of a function at N complex points z i (i = 1,..., N), we define the continued fraction C N (z) = 1 + a 1 a 2 (z z 1 ) a N (z z N 1 ) Here, the coefficients a i are to be determined so that (3.78) C N (z i ) = u i, i = 1,..., N (3.79) The coefficients are then given by the recursion a i = g i (z i ), g 1 (z i ) = u i, i = 1,..., N g p (z) = g p 1(z p 1 ) g p 1 (z), p 2 (3.80) (z z p 1 )g p 1 (z) This requires the following calculation: g 1 (z 1 ) = a 1 = u 1 g 1 (z 2 ) = u 2 g 1 (z 3 ) = u 3 g 1 (z 4 ) = u 4... g 2 (z) = g 1(z 1 ) g 1 (z) g 2 (z 2 ) = a 2 = a 1 g 1 (z 2 ) (z z 1 )g 1 (z) (z 2 z 1 )g 1 (z 2 ) = a 1 u 2 (z 2 z 1 )u 2 g 2 (z 3 ) = a 1 g 1 (z 3 ) (z 3 z 1 )g 1 (z 3 ) = a 1 u 3 g 2 (z 4 ) = a 1 u 4... (z 3 z 1 )u 3 (z 4 z 1 )u 4 g 3 (z) = g 2(z 2 ) g 2 (z) (z z 2 )g 2 (z) g 3 (z 4 ) = a 2 g 2 (z 4 ) (z 4 z 2 )g 2 (z 4 )... g 4 (z) = g 3(z 3 ) g 3 (z) (z z 3 )g 3 (z) g 3 (z 3 ) = a 3 = a 2 g 2 (z 3 ) (z 3 z 2 )g 2 (z 3 ) g 4 (z 4 ) = a 4 = a 3 g 3 (z 4 ) (z 4 z 3 )g 3 (z 4 )... (3.81) Thus, the following triangular matrix p i, j has to be calculated: j = 1 j = 2 j = 3 j = 4... i = 1 a 1 = u 1 u 2 u 3 u 4... i = 2 a 2 g 2 (z 3 ) g 2 (z 4 )... i = 3 a 3 g 3 (z 4 )... i = 4 a
15 This can be done as follows: p 1, j = u j, j = 1,..., N p i, j = p i 1, i 1 p i 1, j (z j z i 1 )p i 1, j, j = 2,..., N and i = 2,..., j (3.82) The diagonal of this matrix then contains the coefficients a i which are needed for the recursion formula for C N (z): with C N (z) = A N(z) B N (z) A n+1 (z) = A n (z) + (z z n ) a n+1 A n 1 (z) B n+1 (z) = B n (z) + (z z n ) a n+1 B n 1 (z) and A 0 = 0, A 1 = a 1, B 0 = B 1 = 1 (3.83) An important fact to keep in mind: the bigger you chose N, i.e. the more coefficients a i you calculate, the bigger the numbers A N (z) and B N (z) are going to be for a given z. N = 128 will already lead to values of order Comparing functions on real and imaginary axes For a first check if an analytic continuation makes sense, there are exact relations at zero frequency. If a complex function is given around z = 0 as f(z) = α + βz + γz (3.84) with complex coefficients α = α + iα etc., we can write it in terms of imaginary frequencies iω n as f(iω n ) = α + iα β ω n + iβ ω n γ ω 2 n iγ ω 2 n +... (3.85) or of real frequencies as f(ω) = α + iα + β ω + iβ ω + γ ω 2 + iγ ω (3.86) and thus by comparing coefficients we find for the function value at zero f(iω n ) ωn =0 = f(ω) ω=0 (3.87) For the derivatives of f we find d Re f(iω n ) dω = d n ωn =0 dω Im f(ω) ω=0 (3.88) 43
16 d Im f(iω n ) dω = d n ωn =0 dω Re f(ω) ω=0 (3.89) d 2 dω 2 Re f(iω n ) = d2 n ωn =0 dω 2 Im f(ω) ω=0 (3.90) d 2 dω 2 Im f(iω n ) = d2 n ωn =0 dω 2 Re f(ω) ω=0 (3.91) Kramers Kronig relations If the function we try to calculate is analytic in the upper complex plane as well as in the lower complex plane (this is the case for all response functions of physical systems) it must obey the Kramers-Kronig relation g (y) = P dx g (x) π y x (3.92) where the notation g(z) = g (z) + ig (z) is used. x, y are used for real, z and u for complex variables. This is the real-axis analog of the more general Cauchy integral formula g(z) = 1 2πi du g(u) u z (3.93) which says that a holomorphic function defined on a disk is completely determined by its values on the boundary of the disk. Response functions are not analytic on the real axis, but otherwise we can deform the contour arbitrarily on the complex plane. If g(z) falls off in infinity, we can write g(z) = which yields g(z) = dx [g(x + iδ) g(x iδ)] 2πi x z dx g (x) π z x The real part of this equation is equal to Eq. (3.92). (3.94) (3.95) The practical calculation of a Kramers Kronig transformation is straightforward but requires the subtraction of the divergent part. It will usually be implemented on a mesh from a to b large enough so that g (a) and g (b) are negligibly small. 44
17 Then g (y) = b a dx g (x) g (y) π x y + g (y) π Re [log(b y) log(a y)] = b a dx g (x) g (y) π x y + g (y) π ( b y ) log y a (3.96) (3.97) Now the integrand is smooth everywhere. At x = y it should be estimated using the derivative dg (y)/dy. Kramers Kronig relations can become very useful if real or imaginary part of a Greens function are easier than the full function. Typically, the imaginary part falls off abruptly at the band edges and is exactly zero beyond. Then the imaginary part is perfectly suited for calculation on a mesh, but the corresponding real part only falls off as 1/ω. 45
where β = 1, H = H µn (3.3) If the Hamiltonian has no explicit time dependence, the Greens functions are homogeneous in time:
3. Greens functions 3.1 Matsubara method In the solid state theory class, Greens functions were introduced as response functions; they can be used to determine the quasiparticle density of states. They
More informationLecture 1: The Equilibrium Green Function Method
Lecture 1: The Equilibrium Green Function Method Mark Jarrell April 27, 2011 Contents 1 Why Green functions? 2 2 Different types of Green functions 4 2.1 Retarded, advanced, time ordered and Matsubara
More information7. Response functions
7. Response functions 7.1 Random Phase Approximation The truncation of the equation of motion for the Greens function with the Hartree Foc mean field approximation transforms the Hamiltonian into an effective
More information2.1 Green Functions in Quantum Mechanics
Chapter 2 Green Functions and Observables 2.1 Green Functions in Quantum Mechanics We will be interested in studying the properties of the ground state of a quantum mechanical many particle system. We
More information561 F 2005 Lecture 14 1
56 F 2005 Lecture 4 56 Fall 2005 Lecture 4 Green s Functions at Finite Temperature: Matsubara Formalism Following Mahan Ch. 3. Recall from T=0 formalism In terms of the exact hamiltonian H the Green s
More informationDrude theory & linear response
DRAFT: run through L A TEX on 9 May 16 at 13:51 Drude theory & linear response 1 Static conductivity According to classical mechanics, the motion of a free electron in a constant E field obeys the Newton
More informationTheory of Many-Particle Systems
Theory of Many-Particle Systems Lecture notes for P654, Cornell University, spring 2005. c Piet Brouwer, 2005. Permission is granted to print and copy these notes, if ept together with the title page and
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechanics Rajdeep Sensarma sensarma@theory.tifr.res.in Quantum Dynamics Lecture #2 Recap of Last Class Schrodinger and Heisenberg Picture Time Evolution operator/ Propagator : Retarded
More informationDiagrammatic Green s Functions Approach to the Bose-Hubbard Model
Diagrammatic Green s Functions Approach to the Bose-Hubbard Model Matthias Ohliger Institut für Theoretische Physik Freie Universität Berlin 22nd of January 2008 Content OVERVIEW CONSIDERED SYSTEM BASIC
More informationOptics and Response Functions
Theory seminar: Electronic and optical properties of graphene Optics and Response Functions Matthias Droth, 04.07.2013 Outline: Light absorption by Dirac fermions Intro: response functions The optics of
More informationQUANTUM FIELD THEORY IN CONDENSED MATTER PHYSICS. Wolfgang Belzig
QUANTUM FIELD THEORY IN CONDENSED MATTER PHYSICS Wolfgang Belzig FORMAL MATTERS Wahlpflichtfach 4h lecture + 2h exercise Lecture: Mon & Thu, 10-12, P603 Tutorial: Mon, 14-16 (P912) or 16-18 (P712) 50%
More informationGreen s functions: calculation methods
M. A. Gusmão IF-UFRGS 1 FIP161 Text 13 Green s functions: calculation methods Equations of motion One of the usual methods for evaluating Green s functions starts by writing an equation of motion for the
More information3.3 Energy absorption and the Green function
142 3. LINEAR RESPONSE THEORY 3.3 Energy absorption and the Green function In this section, we first present a calculation of the energy transferred to the system by the external perturbation H 1 = Âf(t)
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationLandau-Fermi liquid theory
Landau-Fermi liquid theory Shreyas Patankar Chennai Mathematical Institute Abstract We study the basic properties of Landau s theory of a system of interacting fermions (a Fermi liquid). The main feature
More informationF(t) equilibrium under H 0
Physics 17b: Statistical Mechanics Linear Response Theory Useful references are Callen and Greene [1], and Chandler [], chapter 16. Task To calculate the change in a measurement B t) due to the application
More information10. Zwanzig-Mori Formalism
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 205 0. Zwanzig-Mori Formalism Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More information3 Constitutive Relations: Macroscopic Properties of Matter
EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and
More informationSupplementary Information for Experimental signature of programmable quantum annealing
Supplementary Information for Experimental signature of programmable quantum annealing Supplementary Figures 8 7 6 Ti/exp(step T =.5 Ti/step Temperature 5 4 3 2 1 0 0 100 200 300 400 500 Annealing step
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More informationCausality. but that does not mean it is local in time, for = 1. Let us write ɛ(ω) = ɛ 0 [1 + χ e (ω)] in terms of the electric susceptibility.
We have seen that the issue of how ɛ, µ n depend on ω raises questions about causality: Can signals travel faster than c, or even backwards in time? It is very often useful to assume that polarization
More informationBilliard ball model for structure factor in 1-D Heisenberg anti-ferromagnets
Billiard ball model for structure factor in 1-D Heisenberg anti-ferromagnets Shreyas Patankar 1 Chennai Mathematical Institute August 5, 2010 1 Project with Prof. Kedar Damle, TIFR and Girish Sharma, Satyasai
More information2 Canonical quantization
Phys540.nb 7 Canonical quantization.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system?.1.1.lagrangian Lagrangian mechanics is a reformulation of classical mechanics.
More informationFinite Temperature Field Theory
Finite Temperature Field Theory Dietrich Bödeker, Universität Bielefeld 1. Thermodynamics (better: thermo-statics) (a) Imaginary time formalism (b) free energy: scalar particles, resummation i. pedestrian
More information( ) in the interaction picture arises only
Physics 606, Quantum Mechanics, Final Exam NAME 1 Atomic transitions due to time-dependent electric field Consider a hydrogen atom which is in its ground state for t < 0 For t > 0 it is subjected to a
More informationTime dependent perturbation theory 1 D. E. Soper 2 University of Oregon 11 May 2012
Time dependent perturbation theory D. E. Soper University of Oregon May 0 offer here some background for Chapter 5 of J. J. Sakurai, Modern Quantum Mechanics. The problem Let the hamiltonian for a system
More information10. Zwanzig-Mori Formalism
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 0-9-205 0. Zwanzig-Mori Formalism Gerhard Müller University of Rhode Island, gmuller@uri.edu Follow
More information09. Linear Response and Equilibrium Dynamics
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 2015 09. Linear Response and Equilibrium Dynamics Gerhard Müller University of Rhode Island, gmuller@uri.edu
More informationLandau-Fermi liquid theory
Landau- 1 Chennai Mathematical Institute April 25, 2011 1 Final year project under Prof. K. Narayan, CMI Interacting fermion system I Basic properties of metals (heat capacity, susceptibility,...) were
More informationA. Field theoretical methods for the partition function, perturbative expansions
IV. QUANTUM FIELD THEORY, GREEN S FUNCTION, AND PERTURBATION THEORY Markus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus
More informationLecture 20: Effective field theory for the Bose- Hubbard model
Lecture 20: Effective field theory for the Bose- Hubbard model In the previous lecture, we have sketched the expected phase diagram of the Bose-Hubbard model, and introduced a mean-field treatment that
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationThe Feynman Propagator and Cauchy s Theorem
The Feynman Propagator and Cauchy s Theorem Tim Evans 1 (1st November 2018) The aim of these notes is to show how to derive the momentum space form of the Feynman propagator which is (p) = i/(p 2 m 2 +
More information3.5 The random-phase approximation
154 3. LINEAR RESPONSE THEORY 3.5 The random-phase approximation Earlier in this chapter, we have demonstrated that many experimentally observable properties of solids can be expressed in terms of two-particle
More informationClassical field theory 2012 (NS-364B) Feynman propagator
Classical field theory 212 (NS-364B Feynman propagator 1. Introduction States in quantum mechanics in Schrödinger picture evolve as ( Ψt = Û(t,t Ψt, Û(t,t = T exp ı t dt Ĥ(t, (1 t where Û(t,t denotes the
More informationQuantum Physics III (8.06) Spring 2008 Final Exam Solutions
Quantum Physics III (8.6) Spring 8 Final Exam Solutions May 19, 8 1. Short answer questions (35 points) (a) ( points) α 4 mc (b) ( points) µ B B, where µ B = e m (c) (3 points) In the variational ansatz,
More informationQCD at finite Temperature
QCD at finite Temperature I Quantum field theory at finite T François Gelis and CEA/Saclay General outline Lecture I : Quantum field theory at finite T Lecture II : Collective phenomena in the QGP Lecture
More informationPhys460.nb Back to our example. on the same quantum state. i.e., if we have initial condition (5.241) ψ(t = 0) = χ n (t = 0)
Phys46.nb 89 on the same quantum state. i.e., if we have initial condition ψ(t ) χ n (t ) (5.41) then at later time ψ(t) e i ϕ(t) χ n (t) (5.4) This phase ϕ contains two parts ϕ(t) - E n(t) t + ϕ B (t)
More informationGreen's Function in. Condensed Matter Physics. Wang Huaiyu. Alpha Science International Ltd. SCIENCE PRESS 2 Beijing \S7 Oxford, U.K.
Green's Function in Condensed Matter Physics Wang Huaiyu SCIENCE PRESS 2 Beijing \S7 Oxford, U.K. Alpha Science International Ltd. CONTENTS Part I Green's Functions in Mathematical Physics Chapter 1 Time-Independent
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationQuantum measurement theory and micro-macro consistency in nonequilibrium statistical mechanics
Nagoya Winter Workshop on Quantum Information, Measurement, and Quantum Foundations (Nagoya, February 18-23, 2010) Quantum measurement theory and micro-macro consistency in nonequilibrium statistical mechanics
More informationSymmetry of the Dielectric Tensor
Symmetry of the Dielectric Tensor Curtis R. Menyuk June 11, 2010 In this note, I derive the symmetry of the dielectric tensor in two ways. The derivations are taken from Landau and Lifshitz s Statistical
More informationExamination paper for TFY4210 Quantum theory of many-particle systems
Department of Physics Examination paper for TFY4210 Quantum theory of many-particle systems Academic contact during examination: Associate Professor John Ove Fjærestad Phone: 97 94 00 36 Examination date:
More informationEuclidean path integral formalism: from quantum mechanics to quantum field theory
: from quantum mechanics to quantum field theory Dr. Philippe de Forcrand Tutor: Dr. Marco Panero ETH Zürich 30th March, 2009 Introduction Real time Euclidean time Vacuum s expectation values Euclidean
More informationa = ( a σ )( b σ ) = a b + iσ ( a b) mω 2! x + i 1 2! x i 1 2m!ω p, a = mω 2m!ω p Physics 624, Quantum II -- Final Exam
Physics 624, Quantum II -- Final Exam Please show all your work on the separate sheets provided (and be sure to include your name). You are graded on your work on those pages, with partial credit where
More informationFinite-temperature Field Theory
Finite-temperature Field Theory Aleksi Vuorinen CERN Initial Conditions in Heavy Ion Collisions Goa, India, September 2008 Outline Further tools for equilibrium thermodynamics Gauge symmetry Faddeev-Popov
More informationLecture 6: Fluctuation-Dissipation theorem and introduction to systems of interest
Lecture 6: Fluctuation-Dissipation theorem and introduction to systems of interest In the last lecture, we have discussed how one can describe the response of a well-equilibriated macroscopic system to
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationLecture 11: Long-wavelength expansion in the Neel state Energetic terms
Lecture 11: Long-wavelength expansion in the Neel state Energetic terms In the last class we derived the low energy effective Hamiltonian for a Mott insulator. This derivation is an example of the kind
More informationQuantum Grand Canonical Ensemble
Chapter 16 Quantum Grand Canonical Ensemble How do we proceed quantum mechanically? For fermions the wavefunction is antisymmetric. An N particle basis function can be constructed in terms of single-particle
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare ttp://ocw.mit.edu 5.74 Introductory Quantum Mecanics II Spring 9 For information about citing tese materials or our Terms of Use, visit: ttp://ocw.mit.edu/terms. Andrei Tokmakoff, MIT
More informationIntroduction to Green functions and many-body perturbation theory
Introduction to Green functions and many-body perturbation theory Last updated 10 April 2014 Contents 1 Motivation 2 1.1 Connection to experiments............................ 2 1.2 Divergences in the standard
More information4.2 Elastic and inelastic neutron scattering
4.2 ELASTIC AD IELASTIC EUTRO SCATTERIG 73 4.2 Elastic and inelastic neutron scattering If the scattering system is assumed to be in thermal equilibrium at temperature T, the average over initial states
More informationNPTEL
NPTEL Syllabus Nonequilibrium Statistical Mechanics - Video course COURSE OUTLINE Thermal fluctuations, Langevin dynamics, Brownian motion and diffusion, Fokker-Planck equations, linear response theory,
More informationwhere P a is a projector to the eigenspace of A corresponding to a. 4. Time evolution of states is governed by the Schrödinger equation
1 Content of the course Quantum Field Theory by M. Srednicki, Part 1. Combining QM and relativity We are going to keep all axioms of QM: 1. states are vectors (or rather rays) in Hilbert space.. observables
More informationin-medium pair wave functions the Cooper pair wave function the superconducting order parameter anomalous averages of the field operators
(by A. A. Shanenko) in-medium wave functions in-medium pair-wave functions and spatial pair particle correlations momentum condensation and ODLRO (off-diagonal long range order) U(1) symmetry breaking
More informationAttempts at relativistic QM
Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and
More informationIntroduction to Path Integrals
Introduction to Path Integrals Consider ordinary quantum mechanics of a single particle in one space dimension. Let s work in the coordinate space and study the evolution kernel Ut B, x B ; T A, x A )
More informationTopological Work and the Laws of Thermodynamics
Topological Work and the Laws of Thermodynamics Yiheng Xu, 1 Ferdinand Evers, 2 and Charles A. Stafford 1 1 Department of Physics, University of Ariona, 1118 East Fourth Street, Tucson, AZ 85721 2 Institut
More informationEnsembles and incomplete information
p. 1/32 Ensembles and incomplete information So far in this course, we have described quantum systems by states that are normalized vectors in a complex Hilbert space. This works so long as (a) the system
More informationTheoretical aspects of Solid State Physics
Theoretical aspects of Solid State Physics Oleg M. Braun Institute of Physics National Academy of Sciences of Ukraine Contents Introduction 7. The main ideas..................................... 7 2 Mathematical
More informationPhysics 139B Solutions to Homework Set 4 Fall 2009
Physics 139B Solutions to Homework Set 4 Fall 9 1. Liboff, problem 1.16 on page 594 595. Consider an atom whose electrons are L S coupled so that the good quantum numbers are j l s m j and eigenstates
More informationADVANCED TOPICS IN THEORETICAL PHYSICS II Tutorial problem set 2, (20 points in total) Problems are due at Monday,
ADVANCED TOPICS IN THEORETICAL PHYSICS II Tutorial problem set, 15.09.014. (0 points in total) Problems are due at Monday,.09.014. PROBLEM 4 Entropy of coupled oscillators. Consider two coupled simple
More information2 The Density Operator
In this chapter we introduce the density operator, which provides an alternative way to describe the state of a quantum mechanical system. So far we have only dealt with situations where the state of a
More information10 Thermal field theory
0 Thermal field theory 0. Overview Introduction The Green functions we have considered so far were all defined as expectation value of products of fields in a pure state, the vacuum in the absence of real
More information+E v(t) H(t) = v(t) E where v(t) is real and where v 0 for t ±.
. Brick in a Square Well REMEMBER: THIS PROBLEM AND THOSE BELOW SHOULD NOT BE HANDED IN. THEY WILL NOT BE GRADED. THEY ARE INTENDED AS A STUDY GUIDE TO HELP YOU UNDERSTAND TIME DEPENDENT PERTURBATION THEORY
More informationWeek 5-6: Lectures The Charged Scalar Field
Notes for Phys. 610, 2011. These summaries are meant to be informal, and are subject to revision, elaboration and correction. They will be based on material covered in class, but may differ from it by
More information9.1 System in contact with a heat reservoir
Chapter 9 Canonical ensemble 9. System in contact with a heat reservoir We consider a small system A characterized by E, V and N in thermal interaction with a heat reservoir A 2 characterized by E 2, V
More informationList of Comprehensive Exams Topics
List of Comprehensive Exams Topics Mechanics 1. Basic Mechanics Newton s laws and conservation laws, the virial theorem 2. The Lagrangian and Hamiltonian Formalism The Lagrange formalism and the principle
More informationto satisfy the large number approximations, W W sys can be small.
Chapter 12. The canonical ensemble To discuss systems at constant T, we need to embed them with a diathermal wall in a heat bath. Note that only the system and bath need to be large for W tot and W bath
More informationM. A. Gusmão IF-UFRGS
M. A. Gusmão IF-UFRGS - 217 1 FIP164-217/2 Text 9 Mean-field approximation - II Heisenberg Hamiltonian in wave-vector space As we saw in Text 8, the uniform susceptibility does not diverge in the case
More informationQUANTUM MECHANICS. Franz Schwabl. Translated by Ronald Kates. ff Springer
Franz Schwabl QUANTUM MECHANICS Translated by Ronald Kates Second Revised Edition With 122Figures, 16Tables, Numerous Worked Examples, and 126 Problems ff Springer Contents 1. Historical and Experimental
More informationPhys 622 Problems Chapter 5
1 Phys 622 Problems Chapter 5 Problem 1 The correct basis set of perturbation theory Consider the relativistic correction to the electron-nucleus interaction H LS = α L S, also known as the spin-orbit
More informationSummary of Mattuck Chapters 16 and 17
Summary of Mattuck Chapters 16 and 17 Tomas Petersson Växjö university 2008-02-05 1 Phonons form a Many-Body Viewpoint Hamiltonian for coupled Einstein phonons Definition of Einstein phonon propagator
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
More informationMASTER'S THESIS FOR STUD. TECHN. GURO KAHRS ROGNSVÅG. Faculty of Natural Science and Technology Department of Physics
Faculty of Natural Science and Technology Department of Physics MASTER'S THESIS FOR STUD. TECHN. GURO KAHRS ROGNSVÅG Thesis started:.6.7 Thesis submitted: 7.3.7 DISCIPLINE: PHYSICS Norsk tittel: Lineær
More informationSECOND QUANTIZATION. Lecture notes with course Quantum Theory
SECOND QUANTIZATION Lecture notes with course Quantum Theory Dr. P.J.H. Denteneer Fall 2008 2 SECOND QUANTIZATION 1. Introduction and history 3 2. The N-boson system 4 3. The many-boson system 5 4. Identical
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationShigeji Fujita and Salvador V Godoy. Mathematical Physics WILEY- VCH. WILEY-VCH Verlag GmbH & Co. KGaA
Shigeji Fujita and Salvador V Godoy Mathematical Physics WILEY- VCH WILEY-VCH Verlag GmbH & Co. KGaA Contents Preface XIII Table of Contents and Categories XV Constants, Signs, Symbols, and General Remarks
More informationd 3 k In the same non-relativistic normalization x k = exp(ikk),
PHY 396 K. Solutions for homework set #3. Problem 1a: The Hamiltonian 7.1 of a free relativistic particle and hence the evolution operator exp itĥ are functions of the momentum operator ˆp, so they diagonalize
More informationProblem Set No. 1: Quantization of Non-Relativistic Fermi Systems Due Date: September 14, Second Quantization of an Elastic Solid
Physics 56, Fall Semester 5 Professor Eduardo Fradkin Problem Set No. : Quantization of Non-Relativistic Fermi Systems Due Date: September 4, 5 Second Quantization of an Elastic Solid Consider a three-dimensional
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationTheoretical Statistical Physics
Theoretical Statistical Physics Prof. Dr. Christof Wetterich Institute for Theoretical Physics Heidelberg University Last update: March 25, 2014 Script prepared by Robert Lilow, using earlier student's
More informationSelf-consistent Conserving Approximations and Renormalization in Quantum Field Theory at Finite Temperature
Self-consistent Conserving Approximations and Renormalization in Quantum Field Theory at Finite Temperature Hendrik van Hees in collaboration with Jörn Knoll Contents Schwinger-Keldysh real-time formalism
More information3 Symmetry Protected Topological Phase
Physics 3b Lecture 16 Caltech, 05/30/18 3 Symmetry Protected Topological Phase 3.1 Breakdown of noninteracting SPT phases with interaction Building on our previous discussion of the Majorana chain and
More information6. Auxiliary field continuous time quantum Monte Carlo
6. Auxiliary field continuous time quantum Monte Carlo The purpose of the auxiliary field continuous time quantum Monte Carlo method 1 is to calculate the full Greens function of the Anderson impurity
More informationComplex Analysis MATH 6300 Fall 2013 Homework 4
Complex Analysis MATH 6300 Fall 2013 Homework 4 Due Wednesday, December 11 at 5 PM Note that to get full credit on any problem in this class, you must solve the problems in an efficient and elegant manner,
More informationIntroduction to Integrability
Introduction to Integrability Problem Sets ETH Zurich, HS16 Prof. N. Beisert, A. Garus c 2016 Niklas Beisert, ETH Zurich This document as well as its parts is protected by copyright. This work is licensed
More informationFrom Particles to Fields
From Particles to Fields Tien-Tsan Shieh Institute of Mathematics Academic Sinica July 25, 2011 Tien-Tsan Shieh (Institute of MathematicsAcademic Sinica) From Particles to Fields July 25, 2011 1 / 24 Hamiltonian
More information5.4 Given the basis e 1, e 2 write the matrices that represent the unitary transformations corresponding to the following changes of basis:
5 Representations 5.3 Given a three-dimensional Hilbert space, consider the two observables ξ and η that, with respect to the basis 1, 2, 3, arerepresentedby the matrices: ξ ξ 1 0 0 0 ξ 1 0 0 0 ξ 3, ξ
More informationIntroduction to Multiple Scattering Theory
Introduction to Multiple Scattering Theory László Szunyogh Department of Theoretical Physics, Budapest University of Technology and Economics, Budapest, Hungary and Center for Computational Materials Science,
More informationPhysics 607 Exam 2. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2
Physics 607 Exam Please be well-organized, and show all significant steps clearly in all problems. You are graded on your work, so please do not just write down answers with no explanation! Do all your
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationProblem 1(a): Since the Hamiltonian (1) is a function of the particles momentum, the evolution operator has a simple form in momentum space,
PHY 396 K. Solutions for problem set #5. Problem 1a: Since the Hamiltonian 1 is a function of the particles momentum, the evolution operator has a simple form in momentum space, exp iĥt d 3 2π 3 e itω
More information7 Three-level systems
7 Three-level systems In this section, we will extend our treatment of atom-light interactions to situations with more than one atomic energy level, and more than one independent coherent driving field.
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationPHY 396 K. Problem set #5. Due October 9, 2008.
PHY 396 K. Problem set #5. Due October 9, 2008.. First, an exercise in bosonic commutation relations [â α, â β = 0, [â α, â β = 0, [â α, â β = δ αβ. ( (a Calculate the commutators [â αâ β, â γ, [â αâ β,
More informationV. ELECTRON GAS, RANDOM PHASE APPROXIMATION
V. ELECTRON GAS, RANDOM PHASE APPROXIMATION Marus Holzmann LPMMC, Maison de Magistère, Grenoble, and LPTMC, Jussieu, Paris marus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/marus/cours.html (Dated:
More information( r) = 1 Z. e Zr/a 0. + n +1δ n', n+1 ). dt ' e i ( ε n ε i )t'/! a n ( t) = n ψ t = 1 i! e iε n t/! n' x n = Physics 624, Quantum II -- Exam 1
Physics 624, Quantum II -- Exam 1 Please show all your work on the separate sheets provided (and be sure to include your name) You are graded on your work on those pages, with partial credit where it is
More informationLecture 2: Open quantum systems
Phys 769 Selected Topics in Condensed Matter Physics Summer 21 Lecture 2: Open quantum systems Lecturer: Anthony J. Leggett TA: Bill Coish 1. No (micro- or macro-) system is ever truly isolated U = S +
More information